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I.  INTRODUCTION 


Kim  and  Thiele  [1]  have  developed  a  hybrid  AS-MM  technique  (asymptotic- 
moment  method  technique)  to  find  the  induced  currents  on  the  surface  of 
scatterers  in  the  intermediate  frequency  region.  Kaye,  Murthy  and  Thiele  [2], 
[3]  have  significantly  improved  on  that  work  by  eliminating  the  moment  method 
regions  thereby  making  the  method  potentially  applicable  to  large  scatterers 
as  well  as  ones  of  modest  size.  In  the  work  of  Kaye,  Murthy  and  Thiele,  the 
MFIE  was  reduced  to  a  sequence  of  integral  equations  for  both  the  "optics" 
currents  and  the  "correction"  currents,  each  of  which  is  solved  by  iteration. 
On  the  other  hand,  in  Kim  and  Thiele  [1],  the  "correction"  currents  were 
solved  by  the  moment  method.  The  present  technique,  then,  is  a  purely 
iterative  technique  wherein  all  the  integral  equations  are  of  the  same  form 
and  are  similarly  solved  via  iteration.  In  this  report,  the  iterative  method 
is  improved  still  further  by  reducing  the  amount  of  computation  required  as 
discussed  below. 

In  this  report  we  present  further  investigations  on  the  iterative 
technique  [2].  In  [3],  the  optics  currents  on  the.  lit  side  are  computed 
assuming  that  the  shadow-side  current  is  zero.  Hence,  any  two  scatterers 
whose  lit-side  geometry  is  the  same  would  have  identical  "optics"  currents  on 
the  lit-side,  notwithstanding  the  fact  that  the  shadow-side  geometry  is 
different.  The  effect  of  the  shadow-side  geometry  is  taken  into  account  wnile 
computing  the  correction  currents  on  the  lit  side.  This  necessitates  tne 
computation  of  higher-order  correction  currents.  Furthermore,  there  are 
situations  when  the  approximation  of  the  shadow-side  current  oy  zero  current 
is  inappropriate.  As  an  example,  consider  the  situation  of  a  square  cylincer 


being  illuminated  by  a  plane  wave  when  the  angle  of  incidence  is  such  that  one 
of  the  faces  in  the  shaow  region  is  near  grazing.  The  current  on  that  face  is 
quite  appreciable  and  cannot  be  taken  to  be  zero.  We,  therefore,  modify  the 
iterative  technique  by  incorporating  the  effect  of  shadow-side  currents  for  a 
scatterer  with  either  surface  discontinuities  or  a  smooth  surface.  To  obtain 
these  initial  currents,  we  developed  closed  form  expressions  for  currents  on  a 
wedge.  For  smooth  bodies,  Fock  theory  [4]  furnishes  the  initial  estimates  of 
the  shadow-side  current.  Thus,  we  incorporate  the  edge  diffraction  theory  and 
Fock  theory  into  the  ansatz  for  the  iterative  technique  and  obtain  the  hybrid- 
iterative  technique.  We  are  also  able  to  recast  the  integral  equations  for 
correction  currents  in  such  a  way  that  the  technique  consists  in  computing  the 
optics  currents  over  and  again.  This  effects  a  considerable  simplification  of 
the  computer  program.  Furthermore,  incorporating  the  initial  currents  into 
the  iterative  technique  removed  the  necessity  of  computing  the  higher  order 
correction  currents.  Indeed,  frequently  the  zeroth-order  optics  currents 
themselves  are  a  sufficiently  accurate  solution. 

In  Section  II,  we  develop  the  integral  equations  and  discuss  their 
solution  by  iteration.  In  Section  III,  we  present  theory  relevant  to  the 
initial  estimates  of  the  shadow-side  currents  for  both  smooth  bodies  and  edged 
bodies.  In  Section  IV,  we  deal  with  some  special  numerical  considerations. 
Section  V  presents  extensive  computations  illustrating  the  theory  developed  in 
Section  II.  Both  edged  bodies,  exemplified  by  a  square  cylinder  and  smooth 
bodies  exemplified  by  circular  and  elliptic  cylinders  are  considered.  Section 
V  contains  summary  of  Che  technique  and  points  out  further  avenues  of 


n.  GENERAL  THEORY 


A.  MAGNETIC  FIELD  INTEGRAL  EQUATION 

Consider  d  perfectly  conducting  body  illuminated  by  a  plane  wave.  The 
induced  surface  current  density  on  the  body  at  an  observation  point  P  (see 
Fig.  1)  can  be  computed  from  the  MFIE: 

13(A)  =  2nx?i^(A)  +  2nx  L[l3]  (1) 

where  the  operator  L  is  defined  by 


LC3]  =  /  ^(A')  X  ^G(r)  ds'  .  (2) 

I 

The  surface  of  the  body  is  denoted  by  z  and  A  and  A'  are  the  observation  and 
source  point  position  vectors,  respectively,  on  z  and  r  =  |A  -  A‘|.  The  prime 
on  the  gradient  operator  indicates  that  the  differentiation  is  performed  on 
the  source  coordinates.  The  bar  through  the  integral  sign  is  used  to  denote 
the  principle  value  integral  over  z.  is  the  incident  magnetic  field  vector 
and  n  is  the  outward  unit  normal  to  the  surface  at  A.  G(r)  is  tne  free-spaca 
Green's  function  given  by 


3 


for  the  two-dimensional  problem.  is  the  zero-order  Hankel  function 

of  the  second  kind  and  3  is  the  free  space  propagation  constant.  The  time 
dependence  is  taken  to  be  exp(jtut)  and  is  suppressed  throughout. 

B.  OPTICS  CURRENTS 

The  first  step  in  solving  the  MFIE  for  the  induced  surface  current 
density  is  to  divide  the  total  surface  Z  of  the  body  into  Z^  and 

which  represents  the  lit  and  shadowed  regions,  respectively.  The  dividing 
line  between  these  two  regions  is  the  geometrical  optics  shadow  boundary  which 
is  defined  as  the  locus  of  points  satisfying  u  •  n  =  0,  where  u  is  the  unit 
vector  in  the  direction  of  propagation  of  the  incident  field.  Then  3(fJ)  can 
be  expressed  as  follows: 

D(A)  =  <s3^(^)  >  (1  -  5)D^(S)  (5) 

where  5  =  1  if  e  and  5=0  if  R  e  E^.  Eq.  (1)  can  now  be  rewritten  as: 

3(fi)  =  2nx  +  2nx  LC?^]  +  2nx  L[l3^]  .  (6) 

In  the  notation  used  in  Eq.  (6)  and  in  subsequent  equations,  the  surface  o~ 
integration  is  identified  by  the  superscript  on  the  current  density.  We  now 
write  out  Eq.  (6)  explicitly  for  ^  e  e^  and  5  £  z^: 

=  2nx  rl.(^)  +  2nx  -r  2nx  L[u^]  '7; 

u^(^)  =  2nx  2nx  LCJ^j  2nx  L[u^J 
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The  problem  of  solving  for  the  induced  surface  current  density  nas  now 
been  transformed  into  one  of  solving  Eqs.  (7)  and  (3)  for  and  3^.  >ie 
compute  these  currents  in  sequential  fashion  in  the  following  way. 

First,  let  be  an  estimate  of  shadow  current.  For  edged  bodies  li;<e 
the  square  cylinder,  this  current  is  obtained  from  the  wedge  diffraction 
theory.  For  smooth  bodies  like  circular  or  elliptic  cylinders,  this  current 
is  obtained  using  Fock  theory.  A  detailed  discussion  of  determining  31  will 
be  considered  in  Section  III.  Thus, 


.  (9) 

Substituting  Eq.  (9)  in  £q.  (7)  and  noting  that  the  current  on  the  lit  side 
would  now  be  an  approximation  to  the  true  current  and  denoting  this 
approximate  current  by  3 

0  p 


This  "optics  current"  is  a  significant  improvement  over  the  classical  physical 
optics  (?0)  current  since  it  takes  into  account  mutual  interaction  of  current 
on  both  lit-  and  shadow-regions  in  addition  to  the  geometrical  optics  (GC) 
field. 

A  shadow-side  "optics  current",  3^p,  is  now  defined  by 
substituting  3^  =  3  ^  in  (3): 


3,^.^)  =  2nx  3.(3)  -  2nx  L[3j]  ‘  2nx  iQ 
op  1  '  '  op  op 
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In  this  expression,  the  main  contribution  of  the  second  term  on  the  right-nand 
side  will  be  -2nx  to  cancel  the  first  term.  The  3^^  obtained  from  £q. 
(11)  would  be  closer  to  the  true  current  than  our  initial  guess  j|.  Indeed, 
the  only  inaccuracy  incurred  being  that  through  the  approximate  value  taxen 
for  the  lit-side  current. 

Now  that  a  better  approximation  to  the  shadow-side  current  is  available 
than  is  a  new  optics  current  on  the  lit  side  may  be  obtained.  This  lit- 
side  current  may  be  used  to  further  improve  the  shadow-side  current.  Thus, 
higher-order  optics  currents  may  be  obtained  until  the  desired  accuracy  is 
achieved.  The  integral  equations  for  nth  order  optics  currents  may  be  defined 
as  follows 


(12) 

(13) 


^op  0  solution  of  Eq.  (11). 

In  the  iterative  technique  described  in  an  earlier  report,  the  optics 
current  of  zeroth  order  were  improved  upon  by  the  addition  of  correction 
currents.  The  technique  described  here  is  mathematically  equivalent  to  tne 
correction  current  ansatz  of  the  iterative  technique  and  at  tne  same  time  ’s 
much  simpler  from  the  programming  point  of  view. 


Thus,  the  hybrid-iterative  technique  computes  the  currents  on  the  lit  and 
shadow  regions  in  sequential  fashion  starting  with  an  initial  estimate  for  tne 
shadow  current.  All  the  integral  equations  are  of  the  same  form  and  are,  in 
fact,  Fredholm  integral  equations  of  the  second  kind.  Hence,  all  these 
integrals  may  be  solved  by  iteration  as  discussed  in  [2J,  [3j. 


8 


J-.  -r 


ril.  INITIAL  CURRENTS 


m 


An  initial  estimate  of  the  current  on  a  scatterer  with  surface 
discontinuities  like  a  square  cylinder  is  obtained  from  a  knowledge  of 
currents  on  a  wedge.  For  a  scatterer  with  smooth  surface,  Fock  theory  may  oe 
used  to  obtain  the  initial  estimate  of  the  shadow-side  current. 

A.  Shadow-Side  Current  on  a  Wedge: 

Consider  the  wedge  shown  in  Figure  2  illuminated  by  a  TE-plane  wave- 
When  the  angle  of  incidence  41^  is  such  that  the  shadow  boundary  is  not  "close" 
to  face  B,  the  current  on  face  0  is  given  by  [5], 
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Og  - 
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0  <  <  (n-l)Tr  - 


nm 
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(H) 


where , 

k_  (X)  = 


'2  *  '2 
—  e  J  e  dt 

TT 

X 


x„  =  /  23pn  •  (cos  1  +  cos  ~)/sin  i 


i  =  propagation  constant 


p  =  distance  from  the  edge 


n  =  (2  -  ^) 


a  =  interior  wedge  angle. 


P 


When  the  angle  of  incidence  is  such  that  a  shadow  boundary  is  close  to  face  B, 
the  current  is  given  by, 

^  sgn(cot4('^)  iC_  tg 

0  <  p  <  « 

(n-I)iT  -  <  (n-l)iT  . 

where,  B  =  (1  -  1) 


Xg  =  /  0p/2  •  2n  cot  1^1* 

41*  *  )/2n 

Thus,  Equations  (14)  and  (15)  express  the  currents  on  the  shadow  side  of 
a  wedge  in  closed  form.  These  expressions  involve  the  well-known  modified 
Fresnel  functions,  K._  (x)  [6]  and  are  easy  to  compute.  Figure  3  presents 
shadow-side  currents  on  wedge  computed  using  the  above  expressions  as  well  as 
the  exact  solution  obtained  from  the  eigenfunction  solution.  Note  tnat  tne 
agreement  is  excellent. 

These  currents  on  the  wedge  may  readily  be  used  to  obtain  an  initial 
estimate  of  the  shadow  current  on  a  square  cylinder.  Consider  the  geometry  of 
a  square  cylinder  shown  in  Figure  4.  Faces  2  and  3  are  in  the  snadow 
region.  Faces  1  and  2  constitute  a  wedge  illuminated  with 


the  angle  of 


Shadow-side  currents  on  wedi)e. 


DISThNCE  from  EDGE(UAUELENGTHS) 


The  current  on  face  2  may  Oe  approximated  by  the 


incidence  being  „ 

I  ,D 

current  on  a  corresponding  wedge  with  a  =  1/2  and  =  o.  q  .  In  a  similar 

I  '  *  j  0 

fashion,  the  current  on  face  3  also  may  be  approximated. 


B.  Shadow-Side  Current  on  a  Smooth  Body: 

From  Fock's  theory  [4],  the  tangential  components  of  the  magnetic 

field  on  the  surface  of  a  smooth  perfectly  conducting  body  are  equal  to  the 
tangential  components  of  the  external  field  mu  1  tipi i ed  by  a  universal 
function  of  the  reduced  distance  Therefore,  for  TE  polarization 


"tg  '  "Ig  3(5) 
where  ^3 

g(0  =  e-35'^ 


(16) 

(17) 


The  expression  for  the  Airy  function  and  the  contour  C  are  given  in  [4]. 
For  the  definition  of  z,  and  let  us  consider  the  case  where  the  perfectly 
conducting  body  is  an  infinite  circular  cylinder  of  radius  'a'  with  its  axis 
oriented  in  the  z  direction  and  a  plane  wave  is  incident  along  the  negative  y 
axis  with  the  magnetic  field  also  in  the  z  direction,  as  shown  in  Figure  5. 
According  to  Goodrich  [7],  the  distance  from  the  shadow  boundary  originally 
given  in  Fock's  theory  has  been  generalized  to  the  arc  lengtn  on  the  surface 
of  the  cylinder  from  the  shadow  boundary.  Therefore,  the  magnetic  field  on 
the  cylinder  surface  h  1s  given  by 


where  Hq  is  the  magnitude  of  the  incident  plane  wave  and  g(c)  is  defined  tne 
same  way  as  above.  S  is  the  propagation  constant  and  and 
lengths  on  the  cylinder  surface  given  by 

^In  "  ^ 

^2n  =  ^  (^2  ^ 

where 

(p^  =  2tt  -  <p  ,  and 

^2  =  ^  ^  (20) 

as  shown  in  Figure  5  with  0  <  <p  <  2^  .  The  reduced  distances  corresponding  to 
these  path  lengths  are  given  by 

2mr ) 

?2n  “  ^‘*'2  ^  ^ 

Since  g(0  decreases  monotonical ly  away  from  the  shadow  boundary,  after  tne 
wave  creeps  along  the  cylinder  surface  more  than  half  of  the  circle,  tne  field 
decays  to  an  insignificant  value.  Therefore,  only  the  first  term  of  tne 
series  is  retained. 

Once  the  field  on  the  surface  is  found,  one  can  ootain  tne  surface 
Current  density  through  the  following  relation,  namely 

3  =  n  X  TT  '22) 
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where  n  is  the  normal  to  the  surface.  Therefore,  the  total  Fock  current  in 
the  shadow  region  is  given  by 


3^  =  -  H^{e  ^g(?^)  +  ^g(52^}  (23) 

where 

=  (•^)  ^  (2it  -  ) ,  and 

^2  =  M  (24) 

The  computed  results  of  the  Fock  current  for  the  circular  cylinder  of 
radii  equal  to  0.2,  1.2  and  3.2  wavelengths  are  given  in  Figures  6,  7  and  8, 
respectively. 

While  the  expressions  given  here  are  specifically  valid  for  the  circular 
cylinder,  Fock  theory  is  valid  for  all  smooth  bodies  and  for  both  TE  and  TM 
polarizations  and  may  be  used  to  obtain  the  shadow-side  current  in. all  sucn 
cases. 
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circular  cylinder 


Fiijure  8.  Shadow- side  Fock  current  on  a  circular  cylinder 


[V.  NUMERICAL  CONSIDERATIONS 


A.  Surface  Discontinuities 

In  evaluating  the  integrals  during  the  iteration  process,  special 
attention  has  to  be  paid  to  surface  discontinuities  which,  if  not  treated 
correctly,  give  rise  to  numerical  errors.  As  an  example,  consider  the  wedge 
shown  in  Fig.  9.  Here  the  edge  of  the  wedge  is  the  discontinuity.  The 
current  on  face  A  is  given  by 

3^(0  =  2nx  ?l.(i)  +  2^x  /  ;3®(a‘)x  .  (25) 

’  8 

Introducing  the  Greens  function  as  given  in  Eq.  (4)  and  carrying  out  the 
vector  products,  this  becomes: 

J^(i)  =  2H.(i)  +  ii  /  Hj^^(5r)  cos'^di'  .  (26) 

Then,  for  field  points  close  to  the  edge  (z-*-0)  and  source  points  close  to  tne 
edge  (z'+O),  the  distance  r  between  them  tends  to  zero  and  the  Hankel  function 
has  the  following  small  parameter  behavior: 

Hp)(3r).j4,.  (27) 

Under  these  conditions,  the  integrand  contains  an  inverse  distance  function 
which  gives  rise  to  large  error  when  integrated  numerically.  However,  tms 
integral  can  be  evaluated  analytically  as  follows.  From  Fig.  9,  it  can  oe 
seen  that 


22 


Cosip  =  sina 


(28) 


/  2  2 

r  -  ^  I'  +  I  -  Zii'cosa 


Hence,  for  i+0  and  0  >  cosijj  »  0  so  that  for  field  points  close  to  the 
edge  source  points  far  from  the  edge  do  not  contribute  to  the  integral  in 
(28).  The  upper  limit  of  integration  can  be  replaced  by  a  finite  value,  and 
Eq.  (28)  can  be  written  with  the  aid  of  Eqs.  (27)  -(29)  as  follows: 


/(2^0)  =  2h’(Q)  -  J®(ii'>0)  / 


0  l'  +2  -  22 ‘2  COSa 


This  is  a  standard  integral  which  yields  the  following  result: 


/(2-0)  =  2H^(0)  -  J®(2‘*0)  . 


This  result  is  general  and  holds  for  the  optics  and  correction  currents 
alike.  Its  importance  lies  in  the  fact  that  the  edges  of  polygonal  structures 
and  edge-like  protrusions  on  surfaces  do  not  have  to  be  truncated,  as  was  done 
in  [8],  or  dealt  with  using  a  different  approach.  Hence,  the  current  at  the 
edge  can  be  computed  without  the  need  for  resorting  to  a  hybrid  approach. 


B.  CONVEX  SURFACES 

For  surfaces  comprised  of  planar  sections,  the  term  2nx  L[0]  in  the  MFIE 
is  identically  zero  for  observation  and  integration  points  on  the  same  planar 
section.  This,  however,  is  not  true  for  convex  surfaces,  in  which  case  the 
integrals  are  singular  due  to  the  singular  behavior  of  the  Greens  function 
and,  hence,  principle  value  integrals  have  to  be  used. 

The  details  of  the  principle  value  calculation  for  the  three-dimensional 
case  are  given  in  [9],  with  the  result  that 


/  X  G(r)  ds'  =  /  3(S*)  X  ^  G(r)  ds’  -  [4it  -  uj 

T«  r  * 


where  E'  =  2  -  e,  e  being  the  neighborhood  of  the  singularity.  r,  is  the 
absolute  value  of  the  solid  angle  subtended  by  the  deformed  surface  (in 
taking  /)  at  the  field  point  S,  and  has  to  be  determined  from  the  geometry  of 
the  surface.  If  E  is  smooth  at  i?,  then  n  =  2t. 

For  the  two-dimensional  case,  it  will  be  shown  that  the  integral  in  the 
MFIE  remains  finite  despite  the  singularity  in  the  Greens  function.  Consider, 
for  example,  the  case  of  a  cylinder  of  circular  cross  section  as  shown  in  Fig. 
5.  For  this  two-dimensional  scattering  problem,  £q.  (7)  is  given  by: 


^  I  K(S,i?‘  )  J(f)d£' 

I' 


where 


K(^,S- )  =  4^  (n'  .  r' )  H5-)(3r) 


V.  NUMERICAL  EXAMPLES 


In  this  section,  we  illustrate  the  technique  described  in  earlier 
sections.  We  deal  with  scatterers  with  surface  discontinuities  as  well  as 
smooth  surfaces.  Scatterers  with  edges  are  exemplified  by  a  square 

cylinder.  Scatterers  with  smooth  surfaces  are  exemplified  by  circular  and 
elliptic  cylinders.  In  all  cases  discussed,  the  incident  wave  is  a 

Transverse-Electric  (TE)  plane  wave.  Work  on  the  Transverse-Magnetic  (TM) 
case  is  partially  done  and  will  be  the  subject  of  another  report  in  the  near 
future. 

A.  SQUARE  CYLINDER 

Two  examples  are  presented  here.  In  both  cases,  the  size  of  the  cylinder 
is  the  same,  viz.,  w  =  3.7X.  However,  the  angles  of  incidence,  9^.,  are 
different.  The  case,  iji.  »  115°,  has  been  dealt  with  in  an  earlier  report 

[3].  We  study  the  same  case  here  using  our  hybrid-iterative  technique.  It 
was  necessary  in  [3]  to  compute  four  orders  of  correction  currents  to  obtain  a 
sufficiently  accurate  result.  With  the  hybrid  iterative  technique,  zerotn 
order  optics  currents  themselves  give  a  sufficiently  accurate  solution.  We 
also  present  the  first-order  optics  currents  for  comparison.  Figures  10  to  lA 
show  the  results  for  w  =  3.7X  and  <<).  =  115°. 

We  also  consider  the  case  when  =  95°.  Note  that  in  this  case,  Face  2 

of  the  square  cylinder  is  almost,  but  not  quite,  in  the  visible  region.  As 

can  be  seen  from  the  final  result,  the  current  on  this  is  quite  appreciaole 


nyurij  10.  MaynUude  of  current  on  syuare  cylinder.  Order  ^  0. 

w  -  3.7X.  4-.  =  116®.  MM  Moment  Method,  HIM-  Hybrid  iterative 
method. 


and  may  not  be  neglected  initially  as  was  done  in  the  iterative  technique. 
Figures  15-19  show  results  corresponding  to  this  case.  In  this  case  it  was 
necessary  to  compute  up  to  the  second-order  optics  currents. 

8.  SCATTERERS  WITH  SMOOTH  SURFACE 
i )  Ci rcul ar  Cyl i nder : 

In  [2],  Kaye,  Murthy  and  Thiele  used  the  iterative  technique  to  compute 
the  surface  currents  on  a  circular  cylinder  of  radius  3.2A.  These 
computations  include  the  first-order  correction  current  and  are  characterized 
by  a  ripple.  For  the  same  cylinder,  we  computed  induced  currents  using  tne 
hybrid-iterative  technique.  These  results  are  presented  in  Figures  20-22. 
Note  that  there  is  no  ripple  and  that  the  agreement  between  our  results  and 
the  exact  eigenfunction  solution  is  excellent.  Furthermore,  we  needed  to 
compute  only  up  to  the  zeroth-order  optics  currents. 

In  order  to  demonstrate  that  our  technique  worxs  well  for  "small"  bodies, 
we  also  consider  a  0.2A  radius  circular  cylinder.  We  compute  botn 
zeroth-order  and  first-order  optics  currents.  These  results  are  shown  in 
Figures  23-27.  Note  the  excellent  agreement  between  our  results  and  the  exact 
results  obtained  using  the  eigenfunction  solution. 

The  computer  code  we  have  written  is  general  and  may  be  used  to  compute 
the  induced  currents  on  scatterers  with  an  arbitrary  cross-section.  To 
demonstrate  this  generality,  we  have  also  considered  an  elliptic  cylinder 
(Figure  28)  with  semi-major  and  semi-minor  axes  being  1.5x  and  l.OX, 
respectively.  Our  results  are  compared  with  those  of  the  method  of  moments 
(MM).  These  results  are  shown  in  Figure  29-31.  Note  the  excellent  agreement 
between  the  results  obtained  from  the  MM  and  hybrid-iterative  tecnnique. 


D(UAUELENGTHS) 


order 


23,  MdijnitucJe  of  zero-order  optics  current  on 


Figure  2^.  Ftidse  of  zero-order  optics  current  on  a  circular  cylinder  of 
radius  0.2X  with  -  90*^  . 


180.00 


Fiijure  26.  Phase  of  first-order  optics  current  on  a  circular  cylinder  of 
radius  0.2X  with  =  90*^  . 


idilius  O.Za  with 


igurs  23.  Geometry  for  the  elliptic  cylinder  of  5emi-mejor  axis  'V'  an 
semi-minor  axis  "b".  Angle  of  Incicenc  plane  wave  is  ecuai  ~o  u 
degrees. 
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Bi static  scattering  pattern  for  an  elliptic  cylinde 


VI.  CONCLUDING  REMARKS 


In  this  report,  we  present  a  new  technique,  the  hyori d-i terati ve 
technique,  to  compute  the  induced  currents  on  an  arbitrary,  perfectly- 
conducting  scatterer.  This  technique  is  obtained  by  incorporating  the  edge 
diffraction  theory  and  Fock  theory  into  the  ansatz  of  the  iterative  technique 
for  solving  the  magnetic  field  integral  equation  developed  earlier  [2]. 
Furthermore,  the  iterative  scheme  has  been  modified  by  computing  higher-order 
"optics  currents"  instead  of  correction  currents.  We  list  the  following 
advantages. 

i)  Incorporating  the  edge  diffraction  theory  and  Fock  theory  eliminates  the 
need  for  computing  higher-order  correction  currents.  That  is,  speed  of 
convergence  is  improved.  Frequently,  zeroth  order  optics  currents  themselves 
furnish  sufficiently  accurate  solutions.  CPU  time  requirements  are, 
therefore,  less  demanding. 

ii)  In  the  iterative  technqiue,  the  shadow-side  current  is  initially  taken  to 
be  zero.  This  is  inappropriate  in  some  situations.  This  restriction  is 
removed  in  this  hybrid  method.  Hence,  a  wider  class  of  problems  may  now  be 
sol ved. 

iii)  Modifications  in  the  iterative  scheme  introduced  by  replacing  the 
evaluation  of  the  correction  currents  by  the  evaluation  of  higher-order  optics 
currents  simplifies  the  computer  code.  As  a  result,  savings  in  CPU  time  are 
effected. 

It  had  already  been  established  that  our  technique  does  not  give  spurious 
currents  for  frequencies  corresponding  to  interior  resonances  of  a 
scatterer.  Furthermore,  as  the  size  of  the  body  is  increased  our  technique  is 


computationally  more  efficient  than  method  of  moments.  As  an  example,  our 
technique  requires  only  one-third  to  one-quarter  the  CPU  time  required  by  the 
MM  to  compute  the  currents  on  a  3.7X  square  cylinder.  For  larger  scatterers, 
the  advantage  in  CPU  time  is  increasingly  higher. 

Even  though  the  technique  has  been  demonstrated  only  for  2-D,  perfectly- 
conducting  scatterers,  it  may  be  extended  to  more  complicated  structures  like 
coated  scatterers  and  3-0  bodies. 
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